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FURTHER  INVESTIGATIONS  USING  THE  NEUTRAL  CODE 


Introduction 

NEUTRAL  is  a  one  dimensional  code  that  vas  written  to  simulate  ions, 
accelerated  by  an  electric  field,  charge-exchanging  in  a  background 
medium.  This  charge-exchange  mechanism  has  been  proposed  by  Prono  et  al.* 
to  explain  the  premature  shorting  of  the  anode-cathode  gap  in  ion  diodes, 
and  is  also  thought  to  occur  in  the  reflex  switch^-^.  Further  details  can 
be  found  in  reference  (5),  which  we  will  denote  here  as  memol.  This  report 
will  essentially  be  a  continuation  of  memol,  but  we  will  consider  problems 
that  are  closer  to  the  actual  physical  situation  being  modeled.  Ue  will 
also  discuss  and  present  a  parameter  study  that  was  done  with  the  NEUTRAL 
code.  Parameter  studies  are  very  important  in  the  charge-exchange  problem, 
as  it  relates  to  the  reflex  switch,  since  no  measurements  of  relevant 
physical  quantities,  such  as  density  or  electric  fields,  exist  at  the 
present  time. 


Comparison  of  Results 


To  complete  the  final  testing  phase  of  the  NEUTRAL  code,  see  memol  for 
more  details,  we  used  the  energy  dependent  cross  section  a(v)  assumed  by 
Prono  et  al.*.  The  form  of  this  cross  section  was 


ct(w) 


(1) 


with  aQ  =  2.0  x  10-1^  cm2  and  w0  =  10  Kev.  Keeping  the  rest  of  the  input 
variables  the  same  as  in  the  standard  case,  see  again  memol  for  specific 
values,  and  using  the  cross  section  given  in  equation  (1),  NEUTRAL  produced 
figures  1-3.  As  in  memol,  the  dashed  curves  in  figures  2  and  3  are  the 
analytical  results  for  direct  comparison.  However,  unlike  memol,  these 
theoretical  solutions  are  those  found  by  Prono  et  al.^  and  are  given  here 
for  convenience 

Manuscript  approved  September  20,  1990 


1 


mF  „.-BaW 

g(V,X)  =  H(X  -  aV)e 

o 


a-CB(X  -  aW) 
B 


(2) 


G  (X) 
v 


v  F 
o 

aB 


1  -BX 
CB  +  6 


-CBX 


(3) 


G(X) 


where 


X  =  x/X 
V  =  w/w0 

a  =  w0/(qEX) 


C 

B 


-all/ 2 
e 

(i  -  C)_1 


(4) 


Note  that  these  above  equations  are  the  corrected  versions  of  the  Prono  et. 
al.  solutions,  given  by  Creedon^.  As  in  memol,  g  is  the  neutral 
distribution  function  and  Gw  is  the  neutral  energy  flux  which  we  will  not 
use  here.  And  using  the  notation  of  Prono  et  al.*,  G  is  the  net  neutral 
flux.  From  figure  1  we  determined  that  the  ion  flux  F  =  4.3  x  10*^  cm~^ 
sec"*.  Figure  2  shows  the  neutral  flux  profile  as  a  function  of  x.  The 
dashed  curve  is  just  the  plot  of  equation  (4)  with  the  ion  flux  F  mentioned 
previously,  and  the  solid  curve  represents  the  simulation  results.  In 
figure  3  we  have  the  neutral  distribution  at  the  vacuum  plasma  interface 
(see  memol).  Once  again  the  dashed  curve  represents  the  analytical  result, 
and  this  is  given  by  equation  (2)  with  x  =  xmax.  Ve  see  from  these  figures 
that  we  have  very  good  agreement  between  theory  and  simulation.  This  also 
shows  that  the  assumptions  made  by  Prono  et.  al.*  in  their  derivations, 
were  quite  reasonable. 
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At  this  point  one  may  be  wondering  about  the  magnitude  of  the  ion  flux 
F,  and  wether  or  not  this  has  any  physical  basis.  All  we  can  say  is  that 
this  flux  was  arbitrarily  chosen  since  there  has  been  no  measurements,  at 
least  to  our  knowledge,  of  the  neutral  flux.  For  the  test  cases  run  so 
far  it  is  not  necessary  that  the  values  used  as  input  reflect  exactly  the 
real  charge-exchange  problem,  since  we  are  just  looking  for  consistency  in 
our  output.  However,  the  use  of  reasonable  estimates  of  physical 
quantities  is  certainly  desirable.  The  same  holds  true  for  densities  and 
electric  fields,  they  should  not  be  taken  too  seriously  since  we  are  just 
scoping  the  code 
out . 

A  Parameter  Study 

Next  we  used  another  version  of  the  NEUTRAL  code  that  incorporates  a 
spatially  varing  background  neutral  gas  given  by 


For  our  purpose  n0  =  1.0  x  IC1^  cm-^,  which  is  just  the  value  used  in  the 
standard  run  case,  see  memol  for  more  details.  The  quantity  xD  is  a 
measure  of  the  width  of  the  density  profile  and  will  be  varied  in  our 
parameter  study.  This  alternate  version  of  the  NEUTRAL  code  also  contains 
a  charge-exchange  cross  section  that  is  derived  from  experimental  data^. 
The  form  of  this  cross  section  is  similar  to  that  given  by  equation  (1).  A 
simple  parameter  study  was  done  in  which  we  varied  xQ  and  observed  the 
average  neutral  energy  at  the  point  x  =  xmax  =0.1  cm,  which  is  just  the 
width  of  the  charge-exchange  region  in  the  reference  standard  case  (see 
memol).  The  result  of  this  study  is  shown  in  figure  4.  Computer  runs  were 
done  for  x0  =  0.005  cm,  0.05  cm,  0.01  cm,  0-1  cm,  0.5  cm,  and  1.00  cm.  As 
one  can  see  the  neutral  energy  increases  quickly  until  x0  =  0.05  cm,  after 


3 


which  the  energy  only  rises  slowly.  Ve  expect  this  type  of  behavior 
because  when  x0  >>  w0/qE,  where  q  is  the  ion  charge,  E  is  the  magnitude  of 
the  electric  field,  and  w0  is  the  "cut  off  energy"  of  the  charge-exchange 
cross  section  (see  equation  (1)),  spatial  effects  are  unimportant.  It  is 
only  when  xQ  ~  wQ/qE  or  less  that  spatial  effects  become  appreciable.  For 
wQ  =  10  Rev,  q  =  4.8  x  10"-^  esu  and  E  =  1  MV/cm,  we  find  that  x0  -  0.01 
cm.  Thus  one  can  see  from  this  example  that  if  one  knew  the  charge- 
exchange  cross  section  with  reasonable  accuracy,  and  the  electric  field 
inside  the  charge-exchange  region,  then  by  measuring  the  average  neutral 
energy  one  could  estimate  the  thickness  of  the  background  neutral  gas  or 
anode  plasma  (see  memol). 


Time  Rising  Electric  Field 

The  NEUTRAL  code  has  also  been  modified  to  simulate  a  linearly  rising 
electric  field  in  time,  in  the  charge-exchange  region.  The  time  dependent 
electric  field  E(t)  is  assumed  to  have  the  form 


and  t  is  the  time  in  seconds.  The  quantity  dE/dt  is  a  constant  and  is  just 
the  slope  of  the  time  varying  electric  field.  This  modification  to  the 
NEUTRAL  code  is  important  since  the  electric  field  inside  the  anode  plasma 
of  the  reflex  switch  is  not  static,  as  assumed  in  memol.  Thus  we  can 
estimate  the  average  energy  of  a  neutral  as  it  enters  the  vacuum  gap  (see 
memol  again  for  more  details)  as  a  function  of  time.  At  this  point  we 
would  like  to  mention  that  a  time  rising  electric  field  is  only  one  part  of 
the  entire  physical  picture  and  this  will  be  discussed  later  on  in  this 
report.  The  NEUTRAL  code  was  run  with  the  electric  field,  given  by 
equation  (6),  with 
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f  ^7  1  =  7.335681  x  1011  sta---°-—  (7) 

l  dt  J  cm  sec 

and  one  half  this  value,  namely 

f  )  =  3.667841  x  1011  (8) 

l  dt  J  cm  sec 

The  first  value  of  dE/df  was  chosen  such  that  at  the  end  of  the  simulation 
the  electric  field  had  the  final  value  of  3333.3333  statvolt/cm,  which  is 
just  the  electric  field  used  in  the  reference  run  of  memol.  Also  the 
background  neutral  gas  and  the  charge-exchange  cross  section  were  assumed 
to  be  constants  in  these  simulations,  with  the  same  values  as  used  in  the 
formerly  mentioned  reference  run.  This  makes  comparison  with  a  simple 
analytical  model  possible.  We  note  that  this  time  varying  field  version  of 
NEUTRAL  also  uses  an  improved  particle  pusher  algorithm.  Figure  5  is  the 
simulation  using  the  value  of  dE/dt  given  by  equation  (7).  The  code 
results  are  given  by  the  circles  and  the  dashed  straight  line  is  an 
approximate  best  line  fit,  extrapolated  to  zero  energy,  at  time  t  =  0. ,  to 
this  data.  The  solid  line,  with  the  dagger  point  markings,  is  the  plot  of 
the  equation 


=  (qE(t)X) 

(1  )(.-<■« 

max 

(9) 


where  <En>(t)  is  the  average  neutral  energy  as  a  function  of  time  at  the 
vacuum  gas  interface,  or  x  =  xmax.  E(t)  is  of  course  the  time  dependent 
electric  field  given  by  equation  (6).  Note  that  equation  (9)  is  just 
(equation  (15)  of  memol,  but  with  the  static  electric  field  replaced  by 
E(t).  This  substitution,  though  somewhat  questionable,  should  be  valid  if 
the  electric  field  rises  so  slowly  over  time  that  the  system  goes  through 
a  succession  of  equilibrium  states,  i.e.,  the  problem  is  quasirtatic.  As 
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lone  can  see  from  figures  5  and  6,  this  does  not  seem  to  be  a  too  bad  of  an 
assumption.  Figure  6  is  the  same  situation  except  that  dE/dt  is  now  given 
by  equation  (8).  This  case  was  run  to  see  if  a  smaller  dE/dt  would  cause 
the  simulation  results  to  approach  the  expression  given  by  equation  (9). 
As  one  can  see  the  code  predicts  nearly  the  same  energy  at  t  =  2.27  ns  as 
equation  (9),  but  for  later  times  the  accuracy  is  about  the  same  as  in 
figure  5,  roughly  +3  Rev.  To  determine  if  the  difference  in  slopes  as  seen 
in  figures  5  and  6  is  caused  by  code  inaccuracies,  we  plotted  the  average 
neutral  energy  as  a  function  of  electric  field,  as  given  in  figure  7.  The 
code  generated  results  are  given  by  the  circles,  and  as  in  figures  5  and  6, 
the  dashed  line  is  an  eyeballed  best  line  fit,  extrapolated  to  zero.  The 
solid  curve  is  the  analytical  result,  at  x  =  xmax,  given  by  equation  (15) 
in  memol,  see  also  equation  (9).  Ve  also  would  like  to  point  out,  before 
ve  continue,  that  the  neutral  energies  from  the  simulation  were  found  by 
running  the  static  electric  field  case  for  different  field  strengths,  as 
was  done  in  memol.  Referring  again  to  figure  7  we  see  that,  for  a  fixed 
time  step,  the  energies  differ  by  about  4  Kev  at  E  =  3.333  x  1(P 
statvolt/cm,  which  is  1  MV/cm.  These  differences  in  energy  may  be  due  to 
the  time  step  used,  which  is  certainly  the  case  when  E  =  6.666  x  10^ 
statvolt/cm,  or  the  problem  may  be  due  to  the  particle  pusher  in  the 
NEUTRAL  code.  Anyway  we  can  be  reasonably  confident  that  the  difference  in 
slopes  that  are  found  in  both  figures  5  and  6  can  be  attributed  to  the  code 
numerics  and  not  some  physical  process. 


Time  Rising  Current  Density 

As  was  mentioned  earlier,  the  inclusion  of  a  time  varying  electric 
field,  though  very  informative,  is  only  half  of  the  full  time  dependent 
problem.  Usually  an  electric  field  in  a  plasma  is  accompanied  by  a 
current,  this  can  be  seen  from  the  simple  Ohm's  Law,  given  by  (in  c.g.s. 
units) 


3  =  o£ 


(10) 
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where  J  is  the  current  density,  E  is  the  applied  electric  field,  o  is  the 
(scaler)  conductivity,  and  the  arrows  denote  vector  quantities.  We  see 
from  this  simple  argument  that  a  linearly  rising  electric  field  in  time 
would  result  in  a  proportionally  rising  current.  A  similar  reverse 
argument,  again  using  (10)  would  imply  a  time  varying  current  density  is 
related  to  a  time  evolving  electric  field,  given  by  I?  =  where  17  is 
called  the  resistivity  and  h  =  1 /a.  From  the  various  measurements  that 

have  been  done  on  the  reflex  switch  one  certainly  observes  a  time  rising 
current,  see  references  1-4.  Now  assuming  for  the  moment  a  uniform 
current  density  J,  the  current  I  is  related  to  J  by  J  =  I/A,  where  A  is 
some  cross  sectional  area.  So  we  conclude  that  the  full  time  dependent 
problem  must  include  a  time  varying  current  density  in  order  to  be  more 
realistic.  From  a  computational  viewpoint,  if  we  input  a  known  linearly 
rising  current  density  and  electric  field,  this  automatically  determines 
the  resistivity  or  conductivity  by  equation  (10).  This  could  be  very 
helpful  in  trying  to  theoretically  determine  the  value  of  K  in  the 
empirical  formula  for  anomalous  gap  closure  in  the  reflex  switch,  as  given 
by  Creedon  et  al.^»6  This  formula  has  the  form 

1  3 

j  2  .  K  Dak  (11) 


wi  th 

•,  dJ 

J  =  3—  =  a  constant 
at 

where  J  is  the  time  rate  of  change  of  the  current  density  and  ts  is  the 
primary  gap  shorting  time,  with  being  the  primary  gap  spacing.  See 

memol  for  a  diagram  of  the  reflex  switch.  At  this  point  in  time 
preliminary  work  has  been  done  to  model  the  effect  of  a  time  rising  current 
density  and  will  be  the  subject  of  future  reports. 
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Summary  and  Conclusions 


Ve  have  used  the  NEUTRAL  code  to  simulate  the  charge-exchange  neutral 
problem  as  discussed  in  Prono  et.  al.^.  There  was  good  agreement  between 
the  two  results.  A  parameter  study  was  also  run  in  which  we  varied  the 
width  of  the  background  gas  region  and  found  the  corresponding  average 
neutral  energy.  Ve  concluded  that  such  parameter  study  would  be  beneficial 
in  determining  the  width  of  the  anode  plasma  region.  Ve  have  also  modified 
the  NEUTRAL  code  to  handle  a  linearly  time  rising  electric  field.  Ve 
compared  these  results  to  an  analytic  expression  which  was  derived  by 
replacing  the  static  electric  field  with  the  time  varying  one.  A 
comparison  of  the  code  and  the  analytical  expression  showed  limited 
agreement.  It  does  seem  howev<  r  that,  from  the  code  results,  a  linearly 
rising  electric  field  causes  a  similar  increase  in  the  average  neutral 
energy.  It  was  also  noticed  that  by  varying  the  electric  field  in  the 
static  case,  NEUTRAL  predicted  lover  than  predicted  average  neutral 
energies.  This  may  be  caused  by  not  choosing  the  correct  time  step  or  may 
even  be  the  fault  of  the  particle  pusher,  this  should  be  investigated 
further . 


Future  Work 

Ve  are  beginning  to  work  on  the  problem  of  a  time  rising  current 
density,  since  such  a  problem  is  a  more  accurate  representation  of  the 
operation  of  a  reflex  switch.  Such  a  simulation  will  provide  useful  input 
to  future  models  describing  the  shorting  of  the  vacuum  gap. 
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